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Abstract. Spergel and Steinhardt have recently proposed the concept of dark matter 
with strong self-interactions as a means to address numerous discrepancies between 
observations of dark matter halos on subgalactic scales and the predictions of the 
standard collisionless dark matter picture. We review the motivations for this scenario 
and discuss some recent, successful numerical tests. We also discuss the possibility that 
the dark matter interacts strongly with ordinary baryonic matter, as well as with itself. 
We present a new analysis of the experimental constraints and re-evaluate the allowed 
range of cross-section and mass. 

1 Why Interacting Dark Matter Particles? 

Cosmological models with a mixture of roughly 35% collisionless cold dark mat- 
ter such as axions, WIMPs, or any other candidate interacting through the 
weak and gravitational forces only, and 65% vacuum energy or quintessence 
match observations of the cosmic microwave background and large scale struc- 
ture on extra-galactic scales with remarkable accuracy Q^j. However, a number 
of discrepancies have arisen between numerical simulations and observations on 
subgalactic scales. Q We will address these now in turn. 

Halo density profiles. Simulations of collisionless cold dark matter halos 
(CCDM) m|J^,01 consistently show cuspy halo density profiles similar to the 
generalized Navarro-Frenk-White || (hereafter NFW) profile 

r 3 

P( r ) = Po a , , s x 3 _ a (1) 
r a (r + r s ) ,i a 

where r s is a scale radius, which is related to the virial radius r/, of the halo 
via the concentration parameter c = rh/r s . This profile has the feature that it 
diverges as p ~ r~ a as r — > 0. 

The exact value for a found in simulations is a matter of some debate - 
NFW originally found a ~ 1, while more recent, higher resolution simulations 
find somewhat steeper divergences, a ~ 1.5. This is unlikely to be an artifact; if 
anything, resolution and discreteness effects lead to simulations underestimating 
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the steepness of the central density profile. Biases can occur when fitting |T]) to 
simulated halos. For example, the fact that there is a finite number of particles 
per halo make the choice of the halo center somewhat arbitrary. Any departures 
from the "true" center will give the illusion that the density profile is less cuspy 
than it actually is. 

Observationally, gravitational lensing in clusters and the presence of disks in 
galactic halos make it possible to infer their mass distributions. Observations of 
low surface brightness galaxies (LSB) and dwarfs, chosen for their small mass-to- 
light ratios, indicate cores of galactic halos with shallow density profiles pPJTo|]. 
Similarly, there is some evidence that the central regions of cluster halos not 
associated with luminous galaxies appear have flat density profiles (p = r° as 
r^O)0. 

More recently, preliminary results from high quality measurements of rotation 
curves in edge-on LSBs with thin disks pll indicate a best fit of a ~ 0.5, more 
mildly divergent than and inconsistent with CCDM profiles. 

Note that only exquisitely well-measured rotation curves simultaneously con- 
strain both the concentration parameter c and the slope a [ fl3| . In practice, ob- 
servations do not resolve these parameters independently. This means that the 
problem of halo shapes can be rephrased in terms of c; observed halos tend to 
have smaller concentration (c ~ 6 — 8) than simulated ones (c ~ 20). The crucial 
point is that whether one focuses on slope or concentration, simulations and 
observations do not agree. 

Excess of small scale structure. One of the most striking features of 
halos in high resolution CCDM simulations is that they are heavily populated 
with small subhalos or subclumps. The scarcity of observed satellite galaxies |l5| 
could be interpreted as another failure of the CCDM paradigm on small scales. 

Even if there are mechanisms which inhibit star formation in small clumps 
and thereby leave them dark, a large population of clumps may heat the galactic 
disk and even endanger its stability Jl5| , ^6|Jl7[ | . 

This is a delicate issue; there clearly is halo substructure in galaxies (satel- 
lites, dwarfs, high velocity clouds). Any proposed solution must reduce the num- 
ber of satellites but guard against wholesale smoothing of the density field on 
small scales. 

Tully-Fisher relation. Further evidence for less centrally concentrated ha- 
los comes from the failure of CCDM models to account for the zero-point of 
the Tully-Fishcr relation. Mo and Mao | f 8 explain the surprisingly tight rela- 
tion between the luminosity of disk galaxies and their maximum rotation speed 
in terms of the gravitational interactions between galaxy disk and halo. Their 
model requires the presence of a dark halo with concentration c ~ 3, much less 
than c ~ 20 which is found in simulations. Significantly, this pertains to all disk 
galaxies - not just LSBs and dwarfs. 

Bar stability. In high surface brightness spiral galaxies, the persistence of 
bars implies low-density cores . 

It is conceivable that these discrepancies are due to problems with current 
simulations, the quality of present observations or the omission of important 
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astrophysical processes in the models. However, if any one of them persists, they 
may be an indication that the dark matter is not collisionless. 

In this paper we argue that all of the discrepancies listed above may be 
resolved if the dark matter has strong self-interactions. If this proves to be the 
case, the study of fine-scale structure on subgalactic (kiloparsec) scales, a subject 
that has not been given much attention by particle physicists in the past, could 
prove to be a remarkably powerful and precise probe of the properties of dark 
matter. 

2 The Physics of Self-Interacting Dark Matter Halos 

Consider a CCDM halo. The particles in its center have small rms speeds (are 
"cold"). As subclumps sink into the larger halo, their outer, hotter, less tightly 
bound parts are tidally stripped and remain in the outer parts of the halo. If the 
cold centers of the clumps survive, they join the cold cusp of the halo. Because 
of the lack of two-body interactions (the gravitational potential is smooth) there 
is no heat transport into the cold CCDM cusp. 

Virialized CCDM halos are not in thermal equilibrium. When interactions 
are present, however, the second law of thermodynamics requires that heat be 
transferred into the cold center from the outer halo to approach equilibrium. 
Initially, the increased rms speeds spread out the halo cusp into a smoother 
core. At the same time as the inner part is heated, the outer part is cooled, 
the cooler material moves in, is scattered and mixed, so that eventually the 
halo becomes isothermal at late times. Particles will be ejected due to two-body 
interactions, the halo loses heat and gravothermal collapse sets in, moderated 
by ongoing accretion of small subclumps, which continue to heat the halo. 

An important fact to keep in mind is that the heat conduction coefficient does 
not have a simple dependence on the cross section. A traditional analysis a la 
Chapman and Enskog po|j expands about thermal equilibrium which is brought 
about by frequent interactions. In this regime, where the mean free path of the 
particles is much shorter than any characteristic length scale of the system, the 



As we will see below, the interactions discussed here are constructed to be in 
a different regime altogether. Here, the mean free path is of the order of the scale 
of the halo. In this regime, namely the small cross-section limit, heat transport 
increases as the cross-section increases. For this reason a numerical implemen- 
tation of self-interacting dark matter is non-trivial and cannot be achieved by 
solving the standard equations of cosmological hydrodynamics. Also, because of 
the non-monotonic dependence on cross-section, a simple interpolation between 
the case of collisionless dark matter and the hydrodynamic limit is not sufficient. 

The presence of scatterings has another important effect: the particle ve- 
locities become randomized. The random component of particle velocities in a 
collisionless fluid can be highly correlated, leading to non-spherical velocity el- 
lipsoids [El]. Collisions erase these correlations, randomizing the motions, which 



heat conduction decreases with increasing cross-section, n{r) oc 
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leads to more spherical velocity ellipsoids. The resulting lack of a preferred di- 
rection should isotropize the inner parts of halos where scattering is important. 

Spergel and Steinhardt jlj postulated a self-interaction between dark matter 
particles with strength which encompassed the range 

= odd_ = g x 10 _25 _ x x 1Q -23 cm 2 Ge y-l =0.5-6 cm 2 g" 1 (2) 

M 6 w 

many orders of magnitude stronger than for CCDM particles. (The range quoted 
here is narrower than originally proposed in [Q because of additional constraints 
discussed below). A convenient dimensionless quantity measuring the strength 
of the interaction at the solar circle is the optical depth r = s/s so i ar where 
Ssoiar = 28.4 cm 2 /g. Unit optical depth corresponds to a mean free path of order 
the diameter of the solar orbit Askpc = 16kpc at the densities of the solar 
neighborhood (pskpc — 10 5 p C rit — 7 x 10~ 25 g/cm 3 ). One can then estimate the 
size of the optically thick region R opt ~ 8rkpc. 

The range of s originally quoted in pj was chosen to cover a span of optical 
depths straddling unity. On the optically thin side, subhalos are destroyed by 
spallation in which particles are scattered out of the subhalo one-by-one almost 
every time there is an impact with a hotter particle from the surrounding halo. 
On the optically thick side, subhalos are destroyed by evaporation in which the 
subhalo absorbs all the kinetic energy from nearly every incoming hot particle 
of the surrounding halo, heats up and expands. Either mechanism suffices for 
the purpose of destroying dwarf halos in galaxies. However, a problem with 
evaporation is that it is so efficient that galaxy halos in massive clusters probably 
cannot survive, which conflicts with observations. Eq. (||) takes account of this 
problem by shaving off the optically thick regime from the range of s found in 

!■ 

3 Numerical Studies 

To summarize, the qualitative effects of introducing self-interacting dark matter 
(SIDM) in the regime of cross-sections in (^) are: 

• Heat transport leading to modified halo profiles and core collapse on time 
scales larger than the Hubble time. 

• Isotropization of the velocity tensor in dense regions and hence more spherical 
shapes of the inner parts of halos. 

• Limited destruction of halo substructure through spallation. 

Intriguingly, the observations we listed in the introduction suggest precisely this 
sort of behavior. 

This has motivated several authors J2^ , |23| , p4] , ^5|j26| , ^7| to study the effects 
of self-interactions numerically. Apart from simulations which are done in the 
hydrodynamic limit |24|,|25|], all of these find that SIDM halos contain smooth 
cores. However, the longevity of these cores and the time-scale for gravitational 
collapse are still a matter of debate. 
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Fig. 1. Comparison of halo profiles for non-interacting (solid) and self-interacting 
dark matter with cross-sections ooo = 10 -24 cm 2 GeV -1 (dashed) and odd = 
10 -23 cm 2 GeV _1 (dot-dashed). Halos are chosen from a descending range of masses, 
with the four largest halos in our simulation shown in the leftmost column Three short 
rays corresponding to logarithmic slopes 0.5, 1 and 1.5 are shown to guide the eye in 
the inner part of the halo 



In particular, Kochanek and White simulate isolated halos and find fast 
evolution towards gravothermal collapse for higher mass halos f26f . Caveats of 
their simulations include that 1) they begin with steep, high-density profiles 
which are already cuspy before the interactions are turned on, which leads to 
fast evolution because the central density starts from artificially high values, 2) 
they do not include accretion, which ought to moderate gravothermal collapse 
and 3) the cross-section they simulate are high. For a Milky Way type galaxy 
(R ~ 10i?disk — 35 kpc, Mq ~ 7 x 10 11 M Q ) their simulations cover a range 
from s ~ 2.5 — 83g/cm 2 . They observe core collapse for cross-sections above our 
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Fig. 2. Histogrammed distribution of inner halo slopes a for non- interacting (solid) 
and self-interacting dark matter with cross-section odd = 10 -24 cm 2 GeV _1 (dashed) 
and ctbd = 1CF 23 cm 2 GeV _1 (dot-dashed). We only include halos with more than 1000 
particles, corresponding to masses greater than 3.6 x 10 9 Mq 



range (g) . Burkert |23fl has performed simulations comparable to Kochanek and 
White's, and his results suggest significantly longer lived cores. Yoshida et al. 
[[27) compute the halo profile of a cluster from cosmological initial conditions. 
They, too, find large cores which persist until today. 

For our simulations, we modified the freely available GADGET tree N-body 
code Q to include scatterings between dark matter particles using Monte Carlo 
techniques similar to The simulations begin with cosmological initial con- 
ditions in the Lambda model currently favored by the majority of observations 
(J?o = 1 = fi m + — 0.7). We follow the evolution of dark matter, repre- 

sented by 128 3 particles in a box with Ah" 1 Mpc on the side and l/i -1 kpc spatial 
resolution. Periodic boundary conditions ensure that we model the continuing 
accretion of halos correctly. 

We simulate interaction strengths of s ~ 0.06, 0.6, 6 and 60cm 2 /g. Out of 
these the smallest value of s produced results that are barely distinguishable 
from collisionlcss dark matter. The largest cross-section considered significantly 
suppresses the formation of halos, and those that do form are much too large 
and diffuse. But for the two intermediate cases our galactic halos consistently 
developed constant density cores. 

Figure [I] compares the profiles of halos from three of our simulations which 
assume different cross-sections but which are otherwise identical, ft is clear that 
for intermediate cross-sections, within the range suggested by Spergel and Stein- 
hardt, core collapse has not set in. Consistent with our previous discussion of 
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the optically thin regime, we find evidence that the stronger the self-interactions 
are, the smoother the halos become. 

Figure ^ shows the histograms of inner slopes measured for all halos with 
more than 1000 particles in each of the 3 simulations. It is clear that in the 
regime we are exploring, increasing the cross-section smoothes the halo profiles. 

These and other results from our simulations will be discussed in more detail 
in a forthcoming publication p9| |. 

4 Could Dark Matter Consist of an Exotic, Neutral, 
Stable Hadron? 

If self-interactions are responsible for resolving the issues of the halo profile 
and galactic substructure, then the ratio of the dark matter-dark matter cross- 
section cjdd to the mass of the dark matter particle must lie in the range given in 
Eq. (||). Curiously, the ratio for ordinary hadrons, such as neutrons and protons, 
lies in roughly this same range, 0(1O -23 cm 2 / GeV) at low energies. One natural 
possibility is that the dark matter belongs to the hidden sector of a supergravity 
or superstring model with hidden gauge interactions analogous to ordinary gluon 
exchange. A second possibility is that the dark matter is composed of exotic 
hadrons — neutral, stable particles which interact through the strong force with 
ordinary quarks and gluons as well as with themselves. In this case, the ratio of 
scattering cross-section to mass for self-scattering (s) might be comparable to 
the ratio for dark matter-nucleon scattering. We should emphasize that the self- 
interaction proposal does not require that there be strong interactions between 
dark matter and ordinary hadrons; it suffices if the dark matter interacts strongly 
with itself. However, the fact that the cross-section required to solve the halo 
profile problem is so similar to the cross-section for the familiar strong interaction 
makes it impossible to resist considering the possibility. 

If the exotic hadron is neutral and stable and does not radiate photons or 
other massless mesons, then it will act like dissipationless, collisional dark mat- 
ter. Possible examples include the quark-gluino bound state, the SO (in some 
models, the lightest supersymmetric stable particle with odd i?-parity and unit 
baryon number) |3(],|3l| or strangelets |52|. Alternatives are non-hadronic ex- 
otics which might also be arranged to have the desired value of ann/m, such 
as gauge singlet mesons |Q, and Q-balls Of these, the first three naturally 
have values of s = a£>£>/m which lie in the desired range, Eq. Another con- 
dition is that the exotic hadron must not bind to light nuclei. Otherwise, the 
particle would efficiently convert hydrogen and helium to exotic nuclei during 
nucleosynthesis and would be found plentifully. Binding to light nuclei may be 
sufficiently suppressed if the particle has quantum numbers that forbid simple 
pion exchange with ordinary nuclei, as might occur if the particle is an isospin 
singlet. Yet a further constraint is that the cosmic abundance must account for 
the observed matter density. If the particle has similar mass and interactions 
to ordinary hadrons, perhaps a natural explanation can be found for why the 
cosmic energy density in exotic hadrons is comparable that of ordinary hadrons. 
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One may suppose that the notion of exotic hadron is easily discounted. 
Surely the particle is already ruled out by existing cosmological and astrophys- 
ical constraints or by existing accelerator experiments or dark matter detectors 
or searches for exotic nuclei. Surprisingly, the story is not nearly so simple. For 
example, if the particle is strongly interacting and similar to a neutron in most 
characteristics, it evades all currently published constraints known to the au- 
thors. During nucleosynthesis, its kinetic energy density is too small to affect 
the cosmic expansion rate and, if it does not bind to light nuclei, it has no effect 
on ordinary nuclear abundances or evades terrestrial searches for unusual nuclei. 
Only an undetectably small fraction of the particles are captured by stars, such 
as the Sun, so they have no significant effect on stellar evolution, supernovae, etc. 
The scattering cross-section for dark matter on ordinary protons and neutrons 
<iDp would presumably be close to the dark matter-dark matter cross-section, 
(Tdd, so one concern is that halo dark matter collides with ordinary baryons 
in the disk, creating excess high energy cosmic rays and destabilizing the disk. 
We find that these considerations do not rule out the exotic hadron scenario, 
either. In cases where the mass and interaction strength is similar to the more 
abundantly produced neutrons, it is difficult to exclude them with accelerator 
experiments. 

Hence, one is led to consider dark matter searches. However, most dark mat- 
ter searches are aimed at weakly interacting particles (WIMPS), such as axions 
and supcrsymmctric ncutralinos, which can penetrate the atmosphere and deep 
underground unimpeded. Consequently, the detectors are placed at ground level 
or in mines. Strongly interacting particles interact in the upper atmosphere and 
will thcrmalize with the air before they reach ground level. For a 1 GeV particle, 
say, the incoming kinetic energy is 1 keV (corresponding to a virial velocity of 
300 km/s) and a local density of 0.4 GeV/cm 3 . If the particle is strongly inter- 
acting, then the average kinetic energy is an undetectable 30 meV by the time 
the particles reach ground level. 

^From these examples, it is clear that the best way to search for strongly 
interacting particles is in balloons or in space with detectors with little or no 
shielding and low energy threshold. A decade ago, Starkman et al. |35| surveyed 
all current observations and experiments and found that the combination of 
cross-section and mass suggested by Eq. (||) is not constrained for masses near 
1 GeV or so, and no new constraints have been published since - until now. 

Figure 3 represents our new analysis of the constraints on strongly inter- 
acting dark matter |j3l| . This plot assumes a spin-independent cross-section 
for scattering of dark matter from nuclei with atomic mass A that is a da = 
A 2 [m re d(A)/m re d{l)] 2 crDp, as considered by Starkman et al. J35|. It is straight- 
forward to modify the plot for spin-dependent cross-sections or cases with dif- 
ferent A-dependence, as might be expected for some of our candidates. The grey 
strip represents the range of aDo/nri given in Eq. (g), which one might imagine 
is roughly similar to the dark matter-proton cross-section crop/m. The plot has 
three new features compared to earlier, similar exclusion plots: 



Self-Interacting Dark Matter 9 




log 10 m (GeV) 



Fig. 3. Plot of the dark matter-proton scattering cross section od p versus dark matter 
particle mass m showing current experimental limits. The constraints from the XQC, 
IMP7/8 and IMAX experiments are new; the constraints from Pioneer and Skylab ex- 
periments have been re-evaluated. The grey region shows the range of dark matter-dark 
matter cross-section odd in Eq. (2). The strip is displayed to study the possibility that 
dark matter consists of exotic hadrons in which odd /tti and ODp/rn are comparable. 
The square marks the values corresponding to a nucleon, which is similar to what is 
expected for some of the most attractive candidates for dark matter. The new results 
appear to rule out the low mass range, but allow masses greater than 10 5 GeV. 



(1) Some important corrections | p8[ have been applied to constraints from older 
experiments (Pioneer 11 and Skylab) compared to Starkman et al. |3qj . The 
corrections reduce the exclusion region at large cross-section; earlier analyses 
had failed to take proper account of shielding. We also argue that the neutron 
star bound in |55| does not apply in the limit of large cross-sections fsif . The 
re-evaluation opens up a range of the grey-striped region for masses between 10 4 
to 10 17 GeV. 

(2) We publish for the first time results of a complete analysis of the IMP7/8 
cosmic ray silicon detector satellite p3|j36| , |38[ and the IMAX balloon-borne SIMP 
(strongly interaction massive particle) detector |3?],|3^,[}9| . 

(3) We present a provisional exclusion region for the X-ray Quantum Calorimeter 
(XQC) experiment Q], an experiment whose first results are not yet published 
[ [LT| and that was not designed with dark matter in mind. Remarkably, this "inad- 
vertent" dark matter detector provides perhaps the most important constraint, 
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ruling out a large, previously unconstrained region of parameter space spanning 
particles in the 0.3 GeV to 10 5 GeV range that includes the SO and low-mass 
strangelets. (A more sophisticated analysis Q is likely to expand the exclusion 
region somewhat.) The key properties of XQC are that it has impressively sen- 
sitive pixels with a low energy threshold (25 eV for a substantial portion of the 
flight, compared to 160 keV for IMP7) a significant detecting area, minimal at- 
mospheric overburden, and negligible material in front of its 1 sr field-of-view. 
The detector consists of 33 pixels of 1 mm 2 area, each with an x-ray absorber 
(7000 Angstroms of HgTe) deposited on 15 microns of silicon attached to sili- 
con detector structure with an implanted thermometer. Both the HgTe and the 
silicon layers are active in detection; silicon dominates the low-mass region and 
HgTe is more important for excluding the high mass region. Although the de- 
tector flew on a sounding rocket at an altitude of 170 km for only a few minutes, 
the low event rate (probably due to x-ray sources) places strong constraints on 
dark matter candidates. Stronger limits can be obtained with subsequent flights. 

5 Conclusion 

We find that recent simulations |27],^9) based on the self- interaction proposal 
of Spergel and Steinhardt with self-interactions of strength 0.1cm 2 /g < s < 
6cm 2 /g lend strong support for the concept, producing results that fit obser- 
vations significantly better than standard collisionless cold dark matter models. 
The favored dark matter candidates, axions and neutralinos, are effectively col- 
lisionless and, hence, are in some considerable jeopardy. The Spergel-Steinhardt 
proposal has stimulated the interesting possibility that dark matter consists of 
particles that interact through the strong force with ordinary matter. Our re- 
evaluation of constraints leads us to conclude that the exotic hadron possibility 
is now ruled out for a substantial range of masses near 1 GeV and cross-sections 
near 10 _24 cm 2 , eliminating some of the most attractive possibilities. At the 
same time, the re-evaluation has re-opened a region encompassing larger masses 
and cross-sections previously thought to be ruled out. 
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